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ABSTRACT 

Temporal analysis of X-ray binaries and Active Galactic Nuclei have shown that hard X-rays react 
to variation of soft ones after a time delay. The opposite trend, or soft lag, has only been seen in a few 
rare Quasi-periodic Oscillations in X-ray binaries and recently for the AGN, 1H 0707-495, on short 
timescales of ~ 10 3 sees. Here, we report analysis of a XMM-Newton observation of Mrk 1040, which 
reveals that on the dominant variability timescale of ~ 10 4 sees, the source seems to exhibit soft lags. 
If the lags are frequency independent, they could be due to reverberation effects of a relativistically 
blurred reflection component responding to a varying continuum. Alternatively, they could be due 
to Comptonization delays in the case when high energy photons impinge back on the soft photon 
source. Both models can be verified and their parameters tightly constrained, because they will need 
to predict the photon spectrum, the r.m.s variability and time lag as a function of energy. A successful 
application of either model will provide unprecedented information on the radiative process, geometry 
and more importantly the size of the system, which in turn may provide stringent test of strong 
general relativistic effects. 

Subject headings: X-rays: galaxies - Galaxies :Seyfert - Galaxies: individual: Mrk 1040 



1. INTRODUCTION 

Active Galactic Nuclei (AGN) are known to show large 
amplitude variability in X-rays. The varia bility broadly 
correlates with the black hole mass (e.g. IQ'Neill et al.l 
2005), but perhaps a more unifying result is that the 
break in the power spectrum depends on the black hole 
mass and luminosity in a similar way for both AGN 
and Galactic black hole systems, providi ng a fundamen- 
tal p lane of these measured quantities (Mc Hardv et al.l 

l2?)om . 

Study of the variability as a function of energy, es- 
pecially time lags between different energy bands, has 
provided insight into the nature of the variability, the 
radiative processes and geometry of these black hole sys- 
tems. For Galactic black hole systems (e.g. Cyg X-l) 
the lags are found to be proportional to the logarithm 
of the energy bin and more impor tantly they depen d 
on the frequency of the variability (|Nowak et al.l [l999h . 
Similarly, all AGN with well measured time lag spec- 
tra show that the variations in the hard band lag be- 
hind those in the soft band and t he size of the hard lag 
decreases with Fourier frequency (iPapadakis et al.ll2001l; 
Vaughan et al.ll2003t iMcHardy et al.ll2004l lArevalo i et al l 
20061: iMarkowitz et al.l 120071 : lArevalo et al.l 120081) . The 
energy dependent time lags in Mrk 110 were found 
to be proportional to the logarithm of the energy 
(|Dasgupta fc Raol 120061) . These observations demon- 
strate the similarity of the nature of X-ray variability 
in the two different mass scales - X-ray binaries and 
AGN. The frequency dependence of the time lags can- 
not be reconciled with a model where they are produced 
by multiple Compton scattering. Instead, cyl indrical in- 
ward propagation of fluctuations in the disk (jLvubarskiil 
1997), can more naturally explain the dependence. The 
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model requires non-localized origin of the hard and soft 
photons wh ich may be achieved either by a transition 
disk region (IMisrall2000r> or spectral hardening in t he in- 
ner regions (|Kotov et afil200ll lATevalo et al.ll2006Tl . 

These time lags reported are hard i.e. the hard X-rays 
react after a delay to soft ones. For Galactic X-ray bina- 
ries (including Neutron stars) , the opposite trend of a soft 
lag has o nly been reported for kHz Quasi-periodic Os- 
cillations (jVaughan et al.l 119981 : iKaaret et all 11999ft and 
for alternate harmonics of a low fr equency Q PO in the 
black hole system GRS 1915+105 (jCuilll9'99D and seem- 
ingly n ever for the broad ba nd continuum variability. Re- 
cently, iFabian et al.l (|2009|) report a soft lag in the AGN 
1H 0707-495, between 0.3-1.0 and 1.0-4.0 keV bands on 
short timescales of ~ 1000 sees. However, the dominant 
variability on longer timescales exhibits typical hard lags. 
Interestingly, the source spectrum can be modeled with 
a complex reflection component (with both Iron K and 
Iron L emission lines) and the soft lag could be due to 
reve rberation effects be tween the continuum and Iron L 
line (|Fabian et al.1 120091 ) . In general, the temporal anal- 
ysis so far undertaken for AGN and Galactic black holes, 
reveal that for the dominant variability, the hard photons 
lag the soft. 

ASCA observat ion of Mrk 1040 reveals a rather typi- 
cal variable AGN (|Q'Neill et al.ll2lMlTurner et al.l[l999[ ) 
with its spectrum showing I ron line emission and p erhaps 
warm absorber signatures (| Reynolds et al.|[l995l) . Here 
we report analysis of a XMM-Newton observation and 
the surprising result is that in the dominant variability 
timescale of 10 4 sees, the source seems to exhibit soft 
lags. 

2. OBSERVATION, ANALYSIS & RESULTS 

XMM-Newton observed Mrk 1040 on 13 February, 
2009 for the duration of 90.9 ksec (observation ID 
0554990101). We processed and filtered the EPIC-pn 
data in a standard way using the Science Analysis Soft- 
ware (SAS) version 10.0 and using the recent calibration 
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Fig. 1. — The unfolded spectrum of Mrk 1040 fitted by an empiri- 
cal model consisting of absorbed dual-Comptonization components 
and some prominent line and edges (see text for details). The ratio 
of the data to the model shows several systematic features but the 
model is a reasonable first approximation at a few percentage level 
to the overall spectrum. 

files. We examined the full field lightcurve above 10 keV 
for flaring particle background which were identified and 
excluded. We used a continuous exposure of 70.4 ks for 
temporal analysis. For spectral analysis we used the 
same length ofdata as taken in temporal study. We ex- 
tracted source and background spectra from a source- 
centered circular region of radius 36" and an off-source 
region, respectively, and using single and double events 
with FLAG=0. We also extracted lightcurves using the 
above criteria and additionally using quadruple events. 
Response matrices and auxiliary files were generated us- 
ing the RMFGEN and ARFGEN tools. The EPIC-pn 
spectrum was grouped using the SAS task specgroup. 

As is the case for many AGN, the time-averaged X-ray 
spectrum of Mrk 1040 as observed by XMM-Newton is 
complex. There is intrinsic neutral absorption, a soft ex- 
cess feature, multiple Iron line emissions, plus edges and 
absorption features, probably due to a warm absorber. 
The complex Iron line features indicate the possibility 
of a complex blurred reflection component. Moreover, 
as shown below, the source exhibited significant secu- 
lar spectral evolution during the observation and hence 
steady state models may not be appropriate to describe 
the time-averaged spectrum. Our motivation here is not 
to fit the spectrum in detail, but to get a first order ap- 
proximation which describes its overall shape. 

We model the spectrum using an intrinsically absorbed 
dual Comptonization model. One Comptonization com- 
ponent to describe the soft excess and the other for 
the hard X-ray continuum. To this basic model we 
include some prominent line emission and edges. The 
Galactic absorption towards the direction of Mrk 1040 
is ~ 6.64 x 10 20 cm -2 and the data requires an addi- 
tional intrinsic neutral absorption of ~ 2.0 ± 0.2 x 10 21 
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Fig. 2. — Top panel: The 640 sec binned X-ray lightcurves of 
Mrk 1040 in the soft, 0.2-2.0 keV (circles) and hard, 2.0-10.0 keV 
(triangles) energy bands. Bottom panel: The ratio of the hard to 
soft count rates. 

cm~ 2 . The soft excess is described by a Comptonized 
component represented by the XSPEC model "nthcomp" 
(jZdziarski et al.lll996tlZvcki et al.lll999D with spectral in- 
dex Ts — 3.1 ± 0.2 and electron temperature kT e — 
0.34 ± 0.04 keV. A second Comptonization model for 
the hard X-ray emission, is represented by the XSPEC 
convolution model "Simpl" ([Steiner et al J 12009) acting 
on the first Comptonization model "nthcomp" . Thus, 
in this description, the first Comptonization component 
which gives rise to the soft excess, is the origin of the seed 
photons for the second Comptonization component giv- 
ing the hard X-ray continuum. The high energy spectral 
index is T = 1.73 ± 0.02 and the fraction of soft excess 
photons getting scattered is / = 0.05 ± 0.03. The Iron 
line emission feature is modeled as a narrow Gaussian at 
6.36 ± 0.02 keV and a broad one centered at 6.6 ±0.1 
keV with a width of 0.4 ±0.1 keV. Two prominent edges 
at 0.66 ± 0.005 and 0.81 ± 0.01 keV with optical depth 
0.72 ± 0.05 and 0.21 ± 0.03 respectively have been in- 
cluded. For a distance of 51 Mpc, the observed and 
unabsorbed luminosities are 1.2 and 1.9 x 10 43 ergs/s 
respectively. The model is formally unacceptable with a 
X 2 /d.o.f = 281/161. This is evident in Figure [1] where 
the unfolded spectrum and the ratio of the data to model 
is shown. There are systematic variations in the ratio 
plot at several energies which may be due to the presence 
of complex emission/ absorption or blurred reflection fea- 
tures and/or spectral variability during the observation. 
However, as can be seen in Figure [T] the model is a good 
first order approximation of the data at a few percentage 
level. The essence of this spectral analysis is that with 
a high energy spectral index of ~ 1.7, moderate intrinsic 
absorption and soft excess, the spectrum of Mrk 1040 is 
typical of AGNs of its class. 
The top panel of Figure [2] shows the background sub- 
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Fig. 3. — The cross-correlation function for different energy 
bands. The reference energy band is 1.0-2.0 keV. From top to 
bottom: (a) 5.0-10.0 (b) 3.0-5.0 (c) 2.0-3.0 (d) 0.8-1.0 (e) 0.5-0.8 
(f) 0.3-0.5. The solid line is the best fit Gaussian curve and the 
dashed vertical lines mark the centroid. 

tracted X-ray 640 sec binned lightcurves of Mrk 1040 in 
soft, 0.2-2.0 keV (circles) and hard, 2.0-10.0 keV (tri- 
angles) energy bands. There is clear strong correlated 
variability in the two bands on ~ 10 4 sec timescales. 
The bottom panel shows that the ratio of the hard to 
soft count rates (i.e. the hardness ratio) significantly de- 
creases from ~ 0.6 to ~ 0.45 during the course of the 
observation, indicating spectral evolution of the source 
on long (~ 10 5 sec) timescales. It is interesting to note 
that the variation in the hardness ratio is clearly uncorre- 
lated with the variability of the lightcurves. This means 
that the primary variability of the source in ~ 10 4 sec 
timescales is not associated with significant spectral vari- 
ation and hence its origin cannot be due to changes in 
the absorption medium which would effect primarily the 
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Fig. 4. — The r.m.s variability and the time lag as a function 
of energy for Mrk 1040 (left) and Mrk 335 (right column). An 
identical analysis has been undertaken for both the sources. The 
errorbars in the time lags are 1-sigma levels estimated by simula- 
tions. 

soft band. 

Lightcurves binned at 64 sees were generated for seven 
energy bands corresponding to 0.3-0.5, 0.5-0.8, 0.8- 
1.0, 1.0-2.0, 2.0-3.0, 3.0-5.0 and 5.0-10.0 keV. Cross- 
correlation analysis was undertaken between the refer- 
ence band 1.0-2.0 keV and others. The XRONOS task 
"crosscor" was used for the analysis. The functions were 
normalized by dividing by the square root of the product 
of the excess variances of the two series (i.e. with "cross- 
cor" flag "norm = 2" ). The entire lightcurve was taken as 
a single segment and the analysis was undertaken in the 
"slow" mode (i.e. not in fast Fourier transform mode). 
In such a case, the task computes errors on the cross- 
correlation by propagating the measurement errors of 
each time bin. Dividing the light-curve into five or more 
segments, gave negative variance for some segments, ren- 
dering the analysis unfeasible. The cross-correlations in 
the delay time range of —2000 to 2000 sees are plotted 
in Figure [3l The lightcurves are well correlated with the 
cross-correlation function peaking at values > 0.8. In 
this delay range, we fit the functions with a Gaussian 
and obtain the best fit centroid. The best fit Gaussians 
are plotted in the figure (solid lines) and their centroid 
values are marked with a horizontal dashed line which is 
taken as the time-delay between the lightcurves. 

We estimate the significance of the measured time- 
delays through simulations. The intrinsic power spec- 
trum of the source is an approximate power-law with 
index ~ 2. Based on this power spectral shap e and u sing 
the method proposed bv iTimmer Sz Koenigl (|1995t l. we 
simulate 200 pairs of lightcurves having the same mea- 
sured time-delay as the observed pair. Measurement er- 
rors were added to each lightcurve and the pairs were 
then subjected to an identical analysis, where the cross- 
correlation function was fitted by a Gaussian function. 
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The root mean square deviation of the centroids of the 
best fit Gaussian functions, was then taken to be the 
1-sigma error on the time-delay. 

The first column of Figure |4] shows the r.m.s and the 
time lag for Mrk 1040 are plotted as function of energy 
bin. Below 2 keV, the time lag decreases with energy 
which means that the variation in the soft photons oc- 
cur after the corresponding variation in the hard ones i.e. 
there is a soft time lag. For energies > 2 keV, the time 
lag increases with energy as in the regular case of hard 
lags. To validate the analysis technique and to make 
a direct comparison with other sources, the right col- 
umn of Figure U show the r.m.s and time lag for Mrk 
335, computed using an identical analysis of its similar 
length XMM-Newton observation. Although the r.m.s 
varies differently with energy, the time lag for Mrk 335 
increases with energy (i.e. hard lags). This is consistent 
with the more detailed a nalysis, including fre quency de- 
pendence, of this source (|Arevalo et al.ll2008h . 

3. DISCUSSION 

Figure [4] shows that the r.m.s decreases with energy 
for Mrk 1040, while it is nearly constant for Mrk 335. 
This may indicate that the variability in Mrk 1040 is 
due to variations in the absorbing medium. Moreover, if 
the absorbing medium (e.g. warm absorber) reacts to a 
change in the hard X-ray continuum with a time delay, 
this could naturally explain the soft lags observed in the 
source. However, as seen in Figure [JJ the hardness ratio 
variation is on longer timescale and uncorrelated with the 
intensity variation. Thus, such a model cannot explain 
the soft lags observed. 

Frequency dependent time lags in X-ray binaries and 
AGN can be explained in terms of fluctuations prop- 
agat ing from the ou ter regions of the disk to the in- 
ner (|Lvubarskiil ll997[ ). To reconcile soft lags, either the 
waves have to propagate outwards or the hard photons 
have to arise from outer regions, both of which seem 
rather physically unrealistic. 

In the most straight forward interpretation, Comp- 
tonization naturally predicts a time lag between energy 
bands to be ~ R/c \og(E 2 /Ei), where R is the size of the 
region and E2/E1 is the ratio of the energies. While the 
energy dependence is consistent with what is observed 
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in black hole binaries and AGN, the Comptonization in- 
terpretation is often ruled out because it is difficult to 
reconcile with the observed frequency dependence of the 
lags. However, it is important to note that Comptoniza- 
tion lag must exist and should manifest at high enough 
frequencies when the wave propagation lag is small. Thus 
at high frequencies the lag should saturate to the Comp- 
tonization lag values. While, this saturation has not been 
detected, the observed time lags already impose a strin- 
gent upper limit on the size of the Comptonizing region, 
R for AGNs. For A rk 564, the time la g variation of 
~ 50 sec \og{E 2 /Ei) (| Arevalo et al.l 120061 ) requires that 
R < 2 x 10 11 cm or < 2GM/c 2 for at 1O 7 M black hole. 

It is not known whether the lags in Mrk 1040 are fre- 
quency dependent and hence a Comptonization origin 
may still be viable. Soft lags due to Comptonization is 
indeed possible, as invoked to explain the soft lags ob - 
served for kHz QPO in X-ray binaries ()Lee et alj [2001). 
A fluctuation in the electron temperature will lead to 
variations in the hard X-rays after a delay. A fraction of 
these hard X-rays may impinge back on the input photon 
producing region and hence affect the soft photons. This 
will lead to a detectable soft lag. Alternatively, the soft 
lag could be due to reverberation of a complex gravita- 
tionally blurred reflection component to var iations of the 
conti nuum, as proposed for 1H 0707-495 (|Fabian et al.l 
120091 ). It is interesting to note that if like for 1H 0707- 
49, the continuum in Mrk 1040 only dominates in the 1-2 
keV band, the time lag variation with energy observed, 
may be naturally explained. The temporal and spec- 
tral signatures of both these models can be quantified, 
although, especially for the reflection scenario, the com- 
plex and non-intuitive effect of light bending needs to be 
taken into account. Both models can be tested and their 
parameters tightly constrained, because they will need to 
sclf-consistently explain the lag and r.m.s versus energy 
as well as the photon spectrum. Successful application 
of either model will provide rich dividends in terms of 
constraining the radiative processes, geometry and more 
importantly the size of the system and provide oppor- 
tunity to test strong General Relativistic effects. Note 
however, that for both scenarios, the time lag should not 
be frequency dependent. 
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